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Abstract 

An alternative super symmetric SO (10) grand unification model with lopsided fermion mass ma- 
trices is introduced. It generates a large solar-neutrino-mixing angle through the neutrinos' Dirac 
\ mass matrix constrained by the SO(10) group structure, avoiding the fine-tuning required in the 

Major ana mass matrix of right-handed neutrinos. The model fits well the known data on masses 
and mixings of quarks and leptons, and predicts a sizable lepton mixing sin 2 2#i3 ~ 0.074, which 
O ■ is significantly larger than that of the original lopsided model. 
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The discovery of neutrino oscillation has opened up a fascinating window for physics 
beyond the standard model. Experimental data on neutrino mass differences and mixings 
help to constrain various theoretical models of new physics. Assuming three light flavors, the 
lepton Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix is characterized by three 
mixing angles (#12, #13, #23) and three CP-violating phases when neutrinos are Majorana 
fermions. The atmospheric and accelerator neutrino data have determined #23 to a good 
accuracy, and the solar neutrino and reactor neutrino experiments have measured # 12 with 
an even better precision pj. These results have already helped to eliminate a large class 
of neutrino mass matrix models in the literature. The CHOOZ reactor experiment has 
found that sin 2 2#i3, if non-zero, should be smaller than 0.1 0. The next generation of 
neutrino experiments under proposal aims to push the limit to sin 2 29 13 ~ 0.01 0,0], which 
undoubtedly will teach us a great deal about the mechanism of neutrino mass generation. 

If small neutrino masses are assumed to arise from the seesaw mechanism |5[, the first 
thing one learns from present data is that the seesaw scale (the mass of the right-handed 
neutrinos) must be very high (close to 10 15 GeV or so). This strongly suggests that the 
seesaw scale may be connected with one of the leading ideas for new physics beyond the 
standard model, i.e., supersymmetric grand unification theory (GUT) according to which all 
forces and matter unify at very short distances corresponding to energies of order 10 16 GeV. 
Since the GUT models unify the quarks and leptons they build in more constraints and have 
better predictive power which can connect neutrino parameters to the well-determined 
quark parameters. 

The most minimal GUT models that incorporate the seesaw mechanism are based on 
SUSY SO(10) since its 16 dimensional spinor representation contains all fermions of the 
standard model along with the right-handed neutrino needed for this purpose as well as the 
fact that it has B — L as a subgroup whose breaking gives rise to masses to the right-handed 
neutrinos. Depending upon which set of Higgs multiplets is chosen to break the B — L 
subgroup of SO(10) ! there emerge two classes of SO(10) models: one that uses 10ff, 126 
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126 H and 120// [7|,|8|, and the other that uses 10//, 16//, 16 H and 45// |2l HDJ, UJ . While 
most of these models are quite successful in fitting and predicting the known experimental 
masses and mixing angles of leptons and quarks, they predict very different values for the 
poorly-known neutrino mixing angle #13 — majority of models with high- dimensional Higgses 
tend to yield #13 close to the current experimental upper bound and majority of those with 
low-dimensional Higgses generally result in a small #13, hence, a small CP violation in the 
lepton sector. Thus it appears that #13 might be an excellent observable to differentiate the 
two classes of SO(10) models. 

Consider, for example, the SO(10) model with low-dimensional Hig gses and the so-called 
lopsided fermion mass matrices proposed by Albright, Babu, and Barr jiol I2I . The lop- 
sidedness built within the Yukawa couplings between the second and third families generates, 
among other interesting physical consequences, the large atmospheric-neutrino mixing angle 
#23 while keeping V c b in the Cabbibo-Kobayashi-Moskawa (CKM) matrix small. The large 
solar- neutrino mixing angle 6 I 12 , however, is generated less elegantly. It is obtained through 
a fine-tuning which not only requires the 2-3 and 3-2 entries in the Majorana mass matrix 
Mr of the right-handed neutrinos to be of order of parameter e appearing in Dirac mass 
matrices of quarks and leptons, but also requires them to be exactly — e [l^- By varying the 
four parameters in the Mr [l3| , the predication of #13 from this model was found to lie in the 
range of 10~ 5 < sin 2 #13 < 10~ 2 . A narrower range of 0.002 < sin 2 #13 < 0.003 is obtained 
when constraints are imposed on the parameter space. If v e disappearance is observed in 
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the next generation of short baseline reactor experiments the original lopsided model 
would be ruled out. 

Given that the lopsided fermion matrix model is one of the most successful GUT theories 
incorporating all the known experimental facts, two obvious questions arise immediately. 
First, is there a more natural way to realize the large solar- neutrino mixing angle without 
fine tuning? And second, if such an alternative model exists, is # 13 consistently small? 
In this paper, we present a modified lopsided model which uses an alternative mechanism 
to generate the solar-neutrino mixing angle. We assume that the right-handed neutrino 
Majorana mass matrix Mr has a simple diagonal structure, and introduce additional off- 
diagonal couplings in the upper-type- quark and neutrino Dirac mass matrices to generate 
1-2 rotation. We found that all the fermion masses and mixing angles can be fitted well in 
the new model. The mixing angle #13, however, is close to the upper limit from the CHOOZ 
experiment and therefore definitely within the reach of next generation reactor experiments. 

Before we present our model for fermion mass matrices, it is instructive to review some of 
the salient features of the SUSY SO(10) model with lopsided fermion mass matrices 
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Through couplings with a set of Higgs multiplets 10#, 16//, 16# and 45# and the constraint 
from the flavor U(l) x Z 2 x Z 2 symmetry, the fermion mass matrices have the following forms, 

frj \ If] \ 

U = e/3 J Mu , N = J -e \ M v , 

\0 -e/3 1 J \0 e 1 J 

f 7] 5 5'e l 

D = ! S a + r/l\M D , L=l 5 -e | M D 

\5'e^ a + e 1 

M R = I -berj e 2 — e I A R , (1) 
\ arj —e 1 J 

where U, D, L, and N denote up-type-quark, down-type-quark, charged lepton, and neutrino 
Dirac mass matrices, respectively, and Mr is the Majorana mass matrix of the right-handed 
neutrinos. As explained in k|> the various entries in the mass matrices come from differ- 



V 


8 


8'e^ 


5 





(r + e 




-6/3 


1 








C 2 7] 2 


—berj 




—berj 


e 2 


T 


arj 


— e 





ent SO(10) invariants in the superpotential, e.g., r\ from 16il6]10#; e from 16 2 16 3 10#45//, 
5,5' from 1 61 162,3 16^16'//; and a from 16216#16 3 16'h. 

The parameter a is of order one, signaling the lopsidedness between the second and third 
families in D and L. This feature leads to a large left-handed neutrino mixing in the PMNS 
matrix and a small left-handed quark mixing shown in the CKM matrix. The parameter e is 
one order-of-magnitude smaller than a and generates the hierarchy between the second and 
third families. In extending to the first family, 5 and 5' were introduced into the D and L. 
The large solar-neutrino mixing angle is from the left-handed neutrino seesaw mass matrix 
which in turn depends on a very specific structure in Mr. 

Since the lepton mixing matrix is defined as 

Upnms = U\U V , (2) 

the large solar mixing angle can either be generated from W L or U v or a combination of both. 
If there is a non-vanishing 1-2 rotation from U u , it can either be generated from the Dirac 
mass matrix of the left-handed neutrinos or from the Majorana mass matrix of the right- 
handed neutrinos or a combination of both. In the following, we focus on the possibilities 
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in which one of the matrices generates a large solar- neutrino mixing angle, keeping in mind 
though that a general situation might involve a mixture of the extreme cases. In the fermion 
mass model in Eq. (1), the large solar-neutrino mixing is induced mainly by the right-handed 
neutrino mass matrix. 

Thus, an alternative possibility is to produce the large solar-neutrino mixing from the 
charged lepton matrix. In fact, in Ref. [1J|, a model was proposed in which both large solar 
and atmospheric neutrino mixings are generated from the lopsided charged-lepton mass 
matrix. The value of sin 2 26 l i3 is again found to be small, 0.01 or less. 

Here we study yet a third possibility of generating a large size 1-2 rotation in the lepton 
mixing from the neutrinos' Dirac mass matrix N. The easiest way to achieve this might be to 
use a lopsided structure in the 1-2 entries of N. However, this is impossible in group theory 
of SO(10). A large rotation, however, can be generated through 1-3 and 2-3 entries without 
affecting, for example, the quark mass hierarchy between the first and second generations. 
Thus we introduce the following modifications of the up-type quark and neutrino mass 
matrices in Eq. (1), 

/ 1] k — p^ 
V - | ,jj 1 Mi . N = I u | M v 

\K + p LU 1 

M R = I b I A fl , (3) 



The symmetric entries u and k in U and N can be generated from the dimension-5 opera- 
tor 16il6j[16//16#]io, and the antisymmetric p entries in U and N are from dimension-6 
operator 16;16.,[16//16 / //]io45//, where the subscript 10 indicate that the spinor Higgses are 
coupled to 10 of SO (10). Because of the modification, the e entries in D and L now must 
be generated from dimension-6 operator 16 i 16j[16//16 / H ]io45//. We assume as in the past 
that 45// Higgs develops a vacuum expectation value (VEV) in the B — L direction. 16// 
and 16// are the Higgs spinors which break the SO(10) to SU(5) by taking the VEV in the 
singlet direction of SU(5). The second pair of 16^ and 16 H develop VEV in 5 and 5 of 
SU(5), respectively, and therefore the operators involving 1Q H and 16^ contribute to up 
and down sectors as weak doublets, respectively. 

Usually a rotation is connected with the mass spectrum. However, in our case the 1-2 
rotation angle from U will be combined with the 1-2 rotation from D to obtain the Cabibbo 
angle 8 C , and a constraint from the up-type quark spectrum must be avoided. Thus, the first 
two families in the U and N cannot be coupled to each other directly, but can be coupled 
indirectly through the third family. The 1-2 rotations in U and N generated from this way 
are proportional to the ratios 7 = (k — p/3)/u and 7' = (k + p)/u, respectively. 

Taking the approximation rj = 0, the dependence of various mass ratios and CKM ele- 
ments on parameters can be seen roughly from the following approximate expressions (the 
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superscript indicates the relevant quantity is at GUT scale) 



0/0 i 2 ° 
m b /m T ~ 1 - g^pye 

m°/m° ~ , 

m° c /m° t ~ (1 + 7 2 V 2 , 
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1 



l + 7 2 3(l + a 2 ) 



C ^ -=L=(- 7 + fce") 
1 + 7 2 

1 



K° ^ -/==f o, 2 ■ n (7 - fee* + VT + ^i!) , (4) 
+ 7 2 3(a 2 + 1) 

where i# and 9 are defined as tie 161 = 3(5 — aS'e 1 ^)/ (ae) and t R = 2>5\/a 2 + l/(ae). The 
expressions for mass ratios in down-type quark and charged lepton sectors are the same 
as those in the original lopsided model. The expressions for m°/m° and elements in CKM 
matrix are new. These approximations allow us to design strategies to fit various parameters 
to experimental data. 

First, we use the up-type quark and lepton spectra and the parameters in the CKM 
matrix to determine 10 parameters a, e, 5, 5', (ft, lu, 7, 77, Mjj and M D . Our best fit yields a 
and e approximately the same as those in the original lopsided model, and thus the successful 
prediction for the mass ratios m°/m° and m^/m" are kept. The two CKM elements \V® S \ 
and together with the CP violation phase 5c p and the constraint on the product ttfR 
from mass ratio m? e /mP , can fix the ti,, £r, 7 and 9. Then u and rj can be fixed from m° c /m® 
and m°, respectively. The down- type quark mass spectrum come out as predictions. 

To see the dependence of the lepton mixing PMNS matrix on various parameters, we 
construct the Majorana mass matrix of left-handed neutrino from the see-saw mechanism 
m v = -NM^N, 

(if I 'a + (k + p) 2 (k + p)uj t](k — p)/a+ (k + p)\ 
(k + p)u uj 2 uj MIj/A r , (5) 

T}(k - p)/a + (k + p) 1 + (« - pf/a + uj 2 /b J 

which depends on the four unknown parameters, 7', Ar, a and b. With parameter a taking a 
reasonably large value, say, order of 0.001 or larger, the 77 dependent terms can be neglected. 
Then one readily sees that the m v matrix can be diagonalized by a 1-2 rotation of angle 9\ 2 
with tan 9\ 2 = 7', and followed by a 2-3 rotation by angle 9% 3 , with 

, on u 2 v /TTy 2 ^ , > 

tan 2V 23 ~ 1 + {k _ p)2/a + ^ 2/6 _ (1 + y2 )w2 • W 
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The neutrino Majorana masses of the second and the third families are 



m v2 = - 
m u3 = - 



(1 + 7'V 2 + y/l + 7 /2 cj (cot20£ 3 - csc2^ 3 )l Ml / A R , 



1 + 7 ,2 K + V 1 + 7' 2 ^ (cot20£ 3 + csc20£ 3 ) M 2 / A R , (7) 



with m v \ = as the result of the approximation rj = 0. Therefore, the present model 
constrains the neutrino mass spectrum as hierarchial, which means that the parameters in 
the light-neutrino mass matrix, the mass eigenvalues and mixings, do not run significantly 
from GUT to low-energy scales. The mass difference Am 2 12 can be used to fix the right- 
handed neutrino mass scale A R . 

Taking into account rotations from matrices m u and L, we arrive at the elements in the 
PMNS matrix, 

jj ( t R 1 1 \ t R a . 

U e 2 = — ; ; COS# 9 o Sinfloo , 

\^/TT^ 3 v /TTy I vT+^; 23 3 ^ 23 



U u3 = . cos H . 7 1 . sin^oo , 

U e3 = — — = cos^ 3 + [ i 1 ' - — 1 - I ^i.^.. (S) 

The data on the solar-neutrino mixing U e 2, together with the ratio of mass differences, 
^ m2 i2/^ m2 23 = m u2/ ( m u3 ~ m 'u2)i can ^ x l' an d ^23; where the latter depends on a combi- 
nation of a and b. Having fixed 7' and parameters in M R , the atmospheric-neutrino mixing 
and U e3 are obtained as predictions. 
We summarize our input and detailed fits as follows. For CKM matrix elements, we 
take \V US \ = 0.224, \V u b\ = 0.0037, \V c b\ = 0.042, and 5cp = 60° as inputs at electro-weak 
scale. With a running factor of 0.8853 for \V u b\, and \V c b\ taken into account, we have 
I Kill = 0-0033 and \V® b \ = 0.037 at GUT scale. For charged lepton masses and up-type 
quark masses, we take the values at GUT scale corresponding to tan/5 = 10 from Ref. (la ]. 
For neutrino oscillation data, we take the solar-neutrino angle to be # so iar = 32.5° and mass 
square differences as Am 2 12 = 7.9 x 10~ 5 eV 2 and Am 2 23 = 2.4 x 10~ 3 eV 2 . The result for 
the 12 fitted parameters is 

a = 1.83, e = 0.1446, 5 = 0.01, 
8' = 0.014, = 27.9°, r] = 1.02 x 10" 5 , 
to = -0.0466, p = 0.0092, k = 0.0191, 
Mu = 82.2 GeV, M D = 583.5 MeV, A R = 1.85 x 10 13 GeV (9) 

There is a combined constraint on a and b, and thus the right-handed Majorana mass 
spectrum is not well determined. As examples, if a = b, a = —2.039 x 10~ 3 ; and if a = 1, 
b = -1.951 x KT 3 . 

We show the result for the down-type quark masses and right-handed Majorana neutrino 
masses (taking a = b) as follows, 

m° = 1.08 MeV, m° = 25.97 MeV, m° b = 1.242 GeV, 

M 1 = 3.77 x 10 10 GeV, M 2 = 3.77 x 10 10 GeV, M 3 = 1.85 x 10 13 GeV . (10) 
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FIG. 1: The predictions of sin 2 # atm and sin 2 26>i 3 against the mass square difference ratio 
Am 2 23 /Am 2 12 . The region of Am 2 23 /Am 2 12 is obtained from the values of Am 2 23 and Am 2 12 
within their 3a limits. 

The predictions for the mixing angles in the PMNS matrix are, 

sin 2 # atm = 0.49, sin 2 29 13 = 0.074 . (11) 

The result for # atm is particularly interesting: Although the lopsided mass matrix model 
is built to generate a large atmospheric- neutrino mixing angle, the charged lepton mass 
matrix alone produces a 2-3 rotation of 63° instead of 45° because of the constraint from the 
lepton mass spectrum. With an additional rotation # 23 ~ 21° fixed mainly from the ratio 
of mass differences Am 2 12 /Am 2 23 , the nearly maximal atmospheric mixing 44.6° comes out 
as a prediction. If one releases the best-fit value of Am 2 12 and Am 2 23 and only imposes the 
3a constraint as 7.1 x 10" 5 eV 2 < Am 2 12 < 8.9 x 10~ 5 eV 2 and 1.4 x 10~ 3 eV 2 < Am 2 23 < 
3.3 x 10~ 3 eV 2 , one would obtain, as shown in Fig. 1, 0.44 < sin 2 # atm < 0.52 which is well 
within the la limit, and 0.055 < sin 2 2#i 3 < 0.110 which, as a whole region, lies in the scope 
of next generation of reactor experiments. 

Finally, we make some remarks on CP phases in the lepton mixing matrix in our model. 
Since we essentially treat all our parameters as real, the CP-violation in the lepton sector 
is essentially absent. One might wonder why the <fi phase in L and D, which generates 
the CP phase in the CKM matrix, does not give contribution to the CP phases in the 
PMNS matrix. The answer is the lopsided structure of the L matrix. In fact, the unitary 
matrix diagonalizing L^L is nearly real, whereas that diagonalizing D^D is not. There is, of 
course, some trivial CP phases due to specific choices of flavor basis. For example, we find 
5qp INS ~ 7T, and some of the Majorana phases are close to tt/2, all of which are believed to 
be artifacts of the model. 

In summary, we have presented an SUSY SO(10) GUT model for the fermion masses 
and mixings, which is developed from the original lopsided model of Albright, Babu and 
Barr [13]. It contains 13 parameters. After fitting them to experimental data, it yields 
a number of predictions. Whenever the experimental data are available, they work well. 
Most interestingly, the model predicts a sin 2 2# 13 around 0.074, which is significantly larger 
than that from any of previous lopsided models. It can surely be tested through the next 
generation of reactor neutrino experiments. It will also have its characteristic predictions 
for lepton flavor violation, leptogenesis as well as proton decay. These issues are currently 
under investigation. 
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